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The implementation of various methods for the determination of 7 through direct CP violation 
arising in the interference of b — > c and b — > u processes in charged as well as neutral B meson 
decays are considered. We show that the inclusion of D° resulting from D*° — » D° + 7r°(7) say 
via B — > KD*° makes a significant difference in the attainable accuracy for 7. Both exclusive 
and inclusive decays of the B ± and B° to states containing D°/D° followed by both inclusive and 
exclusive decays of the D° are discussed. It is shown that with statistics which might be obtained 
at B factories (5 — 10 x 10 8 B-pairs) a la determination of 7 to ~ ±5° may be possible depending 
on the efficiency of reconstruction, backgrounds and the details of the decay amplitudes involved. 
The role of data from a charm factory as well as effects of D° mixing are discussed. Extraction of 
7 with accuracy that is roughly commensurate with the intrinsic theory error of these methods (i.e. 
around 0.1%), which is an important goal, will require > 10 10 B-pairs, namely a Super-B Factory. 

PACS numbers: 12.15.Hh; 11.30.Er; 13.25.Hw 



I. INTRODUCTION 



The B factories at KEK and SLAC have made remarkable progress in many areas of B-physics, in particular in 
the extraction of Cabibbo Kobayashi Maskawa (CKM) Q] parameters crucial to testing the standard model. Indeed, 
the determination of sin 2f3 via B — > 3 jvbKg in such a way that there is no dependence on theoretical assumptions 
promises to usher in a new era of precision tests of the CKM paradigm|^ y|. 

The determination of the other two unitarity angles, a and 7, without theoretical errors still presents a considerable 
experimental challenge. The corresponding CP violating effects are somewhat harder to observe in channels sensitive 
to a and 7. In addition, effects sensitive to a such as B — ► 7T7t are subject to some electroweak penguin contamination. 

In this paper, we will consider the application of several related methods for determining the angle 7 through the 
decay B~ -> K~D° HHH0- The key idea of these methods is that the decay B~ -> K~ D° and B~ -> K~D° can 
interfere if the D° decays to a common hadronic final state as shown in the generic Feynman diagram Fig. ^ where, 
as an example, of a CP-non-eigenstate (CPNES) case, the mode B~ -> A~[D° -> A+tt"] is shown (here we use D 
to indicate a superposition of D° and D°). Similarly for the case of CP-eigenstate (CPES) [4] an example would be 
D°, D° — ► Ksoj, say. These interferences involves the weak phase 7 together with unknown strong phases from the B 
and D decays. A minimum of two modes that are common to the decays of D°, D° are needed to solve for the weak 
phase 7, along with the two strong phases and the suppressed Bt(B~ — > K~D°) that appears to be very difficult to 
measure experimentally because of serious backgrounds 

There are a number of variations of this process which we will consider here. We will show that an effective method 
for the determination of 7 is to allow the A or the D to assume vector or pseudo-scalar forms. Thus the generic 
reaction has four variations B~ — > A~D°, B~ — > A*~D°, B~ — > A _ D*°, and B~ — > A*~D*°. In principle, this set 
of reactions could be further generalized to higher excited states of D or A but we will not further consider this case 
here. Unless the production of such higher states is surprisingly large, the four reactions listed are likely to be of the 
greatest experimental interest. 

In the cases involving D*° the D*° decays to a D° with tt° or 7 and then the D° decays to the hadronic final 
state common to D° and D . One complication that arises in the case of B~ — > A*~D*° is that there are three 
helicity amplitudes which can, in principle, be separated by considering angular correlations between the decays. 
Here, however, we will just consider the inclusive data, that is without attempting to perform the angular analysis 
such as that considered in 

We will also briefly consider the neutral B decays of the same form such as B° -> A*°D°. Here we will only 
consider the self tagging [3, OH cases, where the flavor of the B is determined by having a charged K in the final 
state, for instance B° — > [A* — * A + 7r°]D°. For non-self-tagging modes such as B° —> A5D one must also take into 
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FIG. 1: Feynman diagrams for B~ — » if - [D° — > A' + 7r~]. Diagrams (a) and (b) are b — > c transitions where D° = D° while 
(c) is the b — > u transition with D° = Z)°. Diagram (a) is color allowed while (b) and (c) are color suppressed; D° — > iC + 7r _ is 
Cabibbo allowed, whereas D° — > i(' + 7r~ is doubly Cabibbo suppressed. 

account the oscillation of the B° in which case the main dependence of the results is on S = j3 — a. + ir = 2/3 + 7 as 
considered in [H [H [H Q [l| . 

The key feature of these decays which allows the determination of 7 is the possibility of interference between D° 
and D° channels. This can be accomplished by the selection of the common hadronic final state of the D° decay. The 
formalism given in the appendix shows how to calculate the total rate for both inclusive and exclusive decays of the 
D°. In this case we refer to a decay which is controlled by a single quantum amplitude as exclusive, for instance decays 
to 2 pseudoscalars or a pseudoscalar and a vector would fall into this category. Decays which are the incoherent sum 
of many amplitudes are inclusive. This would include decays to two vectors where no measurement of the polarization 
is made; multibody decays where the data is integrated over all or part of phase space or states which consist of 
several decay modes of different particle content. In the formalism of the appendix, in such cases we introduce an 
additional parameter R which quantifies the degree of coherence associated with the amplitudes composing the 
inclusive decay. This parameter, along with the average strong phase, must also be determined from the data. 

The same methodology can also be applied to inclusive decays of the parent B. for instance, B~ — > K~D a + X 
or the case B~ — > K*~D*° (summed over polarization). Since this approach makes no assumption concerning the 
structure of the amplitude of such inclusive decays it therefore introduces no potential model dependence into the 
analysis. In contrast, the approach of ll6| assumes a resonance structure but, with this modest assumption concerning 
the amplitude, is able to extract more observables since the phase space dependence is fit in detail; they explicitly 
treat the case of B~ — > AT~7r D°, where the D° decays to a CP eigenstate. 

The decays of the D° can be further categorized according the quark level transition involved, it may either be 
Cabibbo allowed (CA), for instance D° — > K~ir + ; sing ly Cabibbo suppressed (SCS), for instance D° -> K+*K~^ 
or Doubly Cabibbo suppressed (DCS) such as D° -> K+tt'. Note that if D° -> X is DCS then D° -> X is CA 
while if D° — ► X is SCS then so is D° — > X. In addition, some two body decays such as D a — > tt + it~, D° — > K + K~ 
and D° — > Kstt are CP eigenstates (CPES). The bulk of these are of the form K$,l + A; since they are CA and in 
particular about 5% of the branching ratio is to negative CPES (CPES-) of the form Kg + X p^| . 

Section II discusses some of the B and D decays that are of interest. In there we also spell out, in some detail, our 
working assumptions. In section III we discuss the expectation for the size of the CP asymmetry in various modes. 
Section IV illustrates extraction of 7 for some of the combination of data sets. In section V we re- visit the feasibility of 
flavor tagging of D a versus D° using their semi-leptonic decays. Section VI deals with DD mixing effects. In section 
VII we briefly discuss the use of neutral (rather than charged) B-mesons, in conjunction with their self-tagging modes 
and direct CP violation to get 7. Multibody B decays are briefly discussed in section VIII. In section IX we briefly 
mention how these determinations of 7 may play out at a super-B factory and we close with a briefly summary in 
section X. Some details on the formalism are relegated to the appendix. 

II. B AND D DECAY MODES 

In this paper, we will compare the results of various determinations of 7 by carrying out \ 2 scans assuming a certain 
level of input data. 

The actual number of events observed at a given integrated luminosity will therefore depend on phases and often 
even on branching ratios which are unknown. To quantify the statistical error, we will assume that a specific fixed 
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Initial B decay 


Subsequent D u decay 


Number of events 


B- -> K-D" 




25 


B~ -> K-D u 




14 


B~ -> K~D U 


K + TT~ + TVK 


106 


B- -> R-D" 


CPES- 


827 



TABLE I: Initial "core data sample" used where the number of events is assumed to be distributed among the given mode and 
its charge conjugate. The corresponding number of events for the three other initial B~ decays: D*°K~, D°K*~ and D*°K*~ 
are assumed to be the same. 

number of events of each type are seen. This will have the side effect of implying a slight variation in Nb, (the number 
of B mesons) x (acceptance) with the input parameters which will generally lie within a reasonable range. There are 
two advantages of adopting this approach: first of all, the x 2 values determined for a given combination of modes will 
be less dependent on the underlying unknown parameters and secondly, the sample calculations should give an idea 
of the precision which can be reached given the quota of events indicated. 

In order to make contact with actual experimental conditions, we will use a simple model for the acceptance of 
various modes. For each observable particle x in the final state we will assign a relative acceptance R x . In addition 
we will assume there is overall reduction in acceptance of R cut due to the relatively hard acceptance cuts which would 
be required to find the signals suggested in these methods. Following the discussion in0,0], the particular values 
we will use in our numerical estimates are: 

R n ± = 0.95 R K ± = 0.8 R n o = 0.5 i? 7 = 0.5 R v ^ 2l = 0.5 

Rcut — -Z — -r (1) 

6 3 

Perhaps the most striking initial manifestation of CP violation in this system will be the case where B~ — > K~D° 
and the D° subsequently decays to K + n~. As pointed out in [5j the large decay rate B~ — > K~D° coupled with 
the Doubly Cabibbo Suppressed (DCS) decay D° — > K + tt~ gives about the same amplitude as the color-suppressed 
decay B~ — ► K~D° coupled with the Cabibbo Allowed (CA) decay D° — > K + ir~ . The result is that the CP violation 
in this combination is potentially large. In addition, this final state is relatively easy to reconstruct so it will likely 
provide the earliest evidence of CP violation from the interference of b — > u and b — > c transitions. 

In this paper, for simplicity, we will assume a core data sample of 100 D° — > K + ir~ events of this sort which 
are distributed equally among four possible parent decays ( B~ — > T)°K~, B~ — > D*°K~, B~ — > T)°K*~ and 
B~ &*°K*-). In particular we will assume an initial data sample as shown in Table [I] Note that in each case the 
number of events is distributed between the given mode and its charge conjugate. 

In order to estimate what this implies concerning the value of Nb, let us assume that the interference between 
the two channels of B~ — » K~[D Q — > K + ir~] is zero; indeed this should be true on average. In Table IVl discussed 
below, we list the branching ratios which we use for various B and D decays. Using these numbers, the sum of the 
two channels neglecting interference is 26 x 10~ 8 . Since only charged B's are being used here, for the needed 25 K ir 
events, it will require Nr = (3 — 6) x 10 8 , where we have included R cu t = [| — |]; corresponding to 300-600 fb^ 1 . To 
estimate the number of events for the other modes we will combine the above acceptances for each of the final state 
particles. 

We will assume that the product of these two quantities is roughly the same between the different combinations of 
K* and D*° in order to give the same number of events in each case, the other modes tend to have larger branching 
ratios but will likely have smaller acceptances. 

The K + ir~ mode alone given in Table ^ is not sufficient to determine 7, in general at least one additional decay 
mode of D° is required. Another prominent DCS mode which could be used is K* + tt~ 0.0. In principle, this is a 3 
body mode and can interfere to some extent with other quasi 2 body modes; however, the K* is sufficiently narrow 
that taking it as a two body mode should be a good approximation. 

In the CA channel, this mode is about 1.5 times that of D° — > K~ir + and we will assume that this is roughly true 
for the overall B~ — » K~[D° — > K* + tt~]. The K* + has three dominant decay modes, each of which has branching 
ratio ~ |: Kltt + , Kgir + and K + ir°. In the case of Kltt + the Kl will be hard to detect and we will assume that 
the acceptance is negligible. For Ksn + the K$ decays | to 7r + 7r~ giving a final state which is likely to have high 
acceptance. There is, however a background from the CA decay D° — ► tt + [K*~ — > Kstt~] so that some cuts may be 
needed. We will assume that this channel has an acceptance reduced by \ due to the cuts giving a total of \ taking 
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CP Eigeristate 


Br [%] 


Est. Acceptance 


rr t~\ TO/ 1 

eff. Br [%] 




1.14 


0.360 


0.410 


Ks[7] — > 7T + 7T _ 7r U ; 7T + 7r _ 7] 


0.38 


0.231 


0.088 




0.73 


0.648 


0.473 


K s [u) ^3tt] 


1.1 


0.324 


0.356 


-K"s[??' — ► 7r + 7r _ 7; 7r + 7r _ 77] 


0.93 


0.188 


0.175 


i^s[0 — » -K" + -R"~] 


0.47 


0.261 


0.123 


Total CPES- 


4.75 




1.625 


ifs [/2(1270) ^tt+tt-] 


0.22 


0.549 


0.121 


7T + 7r~ 


0.14 


0.648 


0.091 




0.412 


0.512 


0.211 


K S K S 


0.03 


0.567 


0.017 


Total CPES+ 


0.58 


0.55 


0.319 



TABLE II: This table consists of a list of some of the CPES decays of the D°. For each mode, the branching ratio (Br) is 
given as the central value in [l^ | and the acceptance estimated as discussed in the text. The effective branching ratio is the 
product of these two. 



into account Br(Ks — > 7r + 7r~). There is no such problem with K* + — » K + ir° but the neutral pion in the final state 
will reduce the acceptance by 0.5. 

Overall then, the acceptance for B~ — > K~[D° — > K* + tt~] is 0.22. The total number of events then should be .55 
times the case of D° — > K + tt~ so we will assume 14 events for each combination as shown in Tabled 

Another class of modes are CP eigenstates (CPES) as considered in £| which we will refer to as GLW modes. Such 
modes, excluding those with Kl have a branching ratio of about 5%. 

In Tabic [H] we list some such modes with the central value of the branching ratio from |18| and acceptance for 
B~ — > K~[D° — > CPES] estimated as above. The product of these two gives an effective branching ratio; for CP=-1 
(CPES-) states it is 1.6% while for CP=+1 (CPES+) this is .32%. In our estimates we will consider the CPES- states 
only in which case if there are 25 K + tt~ events there should be about 827 CPES- events as indicated in Table [i] 

Singly Cabibbo suppressed modes which are not CP eigenstates such as D° — > K* + K~ have been considered in [l7| . 
These have the property that the decay of D° to the charge conjugate i.e. D° — > K + K*~ will also be possible. We 
will therefore be able to observe four different reactions: 



B~ - 


-> K~ 


[D°- 




K~ 


B+ - 


- K+ 


[D°- 


-> K*- 


K+ 


B~ - 


-> K~ 


[D°- 


-> K*- 


K+ 


B+ - 


-+ K+ 


[D°- 


■* K*^ 





Note that the strong phase difference in the first two decays is the opposite of that in the second two. In fact, both of 
these modes taken together will allow the determination of 7. The central values of the two D branching ratios 0] 
are 0.38% for D° — > K* + K~ and 0.22% for D° — > K*~K + . Estimating the number of such events in our data sample 
as above we find about 50 in each case. 

Another data set which might be obtained is by using semi-leptonic (SL) decays of the D° to directly monitor the 
b — > u transition B~ — > if(*)~l)(*) . This is likely to be difficult as discussed in p| as it requires flavor tagging of 
D° versus D° which has non-trivial backgrounds from semi-leptonic decays of the B. In Section V, we will discuss 
some possible strategies for overcoming this potential difficulty. To illustrate possible usefulness of this determination, 
we will then show the effect on 7 extraction assuming that data of this sort gives a 1-a error on Br(B~ — > K~D Q ) 
ranging from 200% to 10%. Thus the implementation of the original GLW0 method to obtain 7, may be envisioned 
with the use of CPES along with flavor tagging through the use of SL data. 

Additional information to reconstruct 7 may also be obtained through the use of a "0(3770) charm factory to measure 
the strong phases in D° decays 0, H3, l2ll |22| . As discussed in this can provide additional information so that 7 
can also be determined through inclusive modes. Even if charm factory data is not available, if we combine inclusive 
modes with other exclusive modes, there is enough information to determine 7 and we shall explore a number of such 
possibilities below. 

In particular, as an example of such a mode we will consider the combination of modes K + n~ + mr. Most likely 
each of these modes will be detected separately so additional information |(| [2j| may be obtained by separating such 
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CP Eigenstate 


Br [%] (CA) 


Est. Acceptance 


eff. Br [%] 




11.1 


.304 


3.37 


K~ 1T + TT + TT~ 


7.5 


.549 


4.12 


K~ tt + it + it~ 7r u 


4.0 


.275 


1.10 


Total 


22.6 




8.59 



TABLE III: The Cabibbo allowed (CA) branching ratio for various decays of the form D° — > K~ir + + nir are shown together 
with the acceptance as discussed above and the effective branching ratio. We will assume that the corresponding doubly- 
Cabibbo-suppressed (DCS) modes are equal to these multiplied by a factor of sin 4 (9c • In the case of K~ir°ir + we have 
subtracted the portion which is produced via the K*~ir + channel to avoid double counting with that mode. 

events according to particle content or by analyzing the phase space distributions. In the early stages, however, such 
analysis may not give significant gains over treating these modes as inclusive ones. 

In Table IIHI we show the acceptance and effective branching ratios for various modes of this form. Note that we 
subtract ir + [K*~ — > K~tt°] from the sample to avoid double counting. Again we will assume that the DCS modes 
are related by a factor of sin 4 6c- Because this is an inclusive state, we also need to specify the coherence factor Rp 
defined in the appendix. We will assume Rp = 0.51 for our numerical illustration which is the value determined in the 
model discussed in Q. The branching ratio to these modes via the CA channel are about 7 times that of D° — > K + ir~ 
with a similar ratio being true for the DCS analogs. Assuming the acceptance is roughly 1/2 of the two body decay, 
the number of events in this sample (106) is thus as shown in Tabled 

In summary, we will investigate the determination of 7 through considering various combinations of the following 
types of data shown in Table [I] 

1. D° -> K+7T-. 

2. D° -> K*+ir-. 

3. GLW modes, i.e. D° -> CPES. 

4. Inclusive decays of the form D° — > K + ir~ + nir. 

5. Singly Cabibbo Suppressed D decays such as D° — > K* ± K^ . 

6. A determination of B~ — > K~D° via a scmi-lcptonic tag of D°. 

7. A determination of strong phases and the coherence factor, Rp, from a "0(3770) charm factory 
These will be combined with the following decays of the initial B meson: 

1. B~ -» if(*)-D(*)° 

2. B~ -> K-T>° +X 

3. B~ -> [K*- -► A—7r°]D + X 

The accuracy in determining 7 will, in general, depend on the value of 7, the strong phase of the B decay and the 
strong phase of the D decays. Clearly, these quantities are unknown and cannot be calculated. We are thus driven 
to consider sample calculations with arbitrarily chosen strong phases in order to project how the measurement will 
proceed. 

For our sample calculations we will assume that 7 = 60° which is consistent with current data (24|. The strong 
phase differences for the various B and D decays we will use in our illustrative calculations are chosen completely 
arbitrarily and are given in Table IIVI The table gives the strong phase difference between the B decay involving 
a D° in the first column and the corresponding decay involving D° in the second column. Note that in the case 
of B~ — > K*~T)*° there are three different helicity amplitudes, helicity h — +1, 0, and — 1, hence there are three 
different strong phases. 

In order to carry forward our analysis we will need to make some assumptions concerning branching ratios of decays 
which have not been measured yet. In Tablc^wc summarize the branching ratios for the decays which we use. The 
decays which are measured are taken from [l8j and indicated by (*) and the ones which we estimate are indicated by 
(t). 
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Decay 1 


Decay 2 


Phase Difference 


B~ -> D"K- 


B~ -> D U K~ 


-50° 


B — > D K 


B~ — * D"*K~ 


— 10° 


i5 — > U ri 


_D — > U J\ 


+oU 




r- _^ n*° if*- 
R~ — > D*° ft"" - 
_o u riL h _ +1 


-1-70° 

+110° 

-4-1 ^0° 


B — > L> A 


B — ► A> A 


130° 


B o ^ D *o^_ i 


B- _» D*°K* h ° =+1 
B- - ZT^" +1 


+50° 
+10° 
-30° 


B" -» B U A~ + A 


B~ D U A~ + X 


45° R = 0.7 


B ,J -» D"K~ + X 


B u -> D u K~ + X 


135° R = 0.7 


D u - A+^- 


D 11 -» A+^~ 


+120° 


D U -> A + TT" + 72 7T 


B U -> A+7T~ + T17T 


+90°(B = 0.5) 




£>" _ A*+ T " 


+60° 


D u -> K*+K~ 




+30° 



TABLE IV: The strong phase differences used in our sample calculations. The strong phase difference is between decay 1 
involving D° and decay 2 involving D°. 



Decay 


Branching Ratio 


Status 


B u - K+7T- 


3.80% 


* 




1.48 x 10~ 4 


* 


B u -» A*+^" 


6.0% 


* 


D u - A*+^" 


1.7 x 10" 4 


t 


B U -> A + 7T~ + 717T 


25% 


t 


B" -> A+7T~ + TOT 


7 x 10" 4 


t 


D u -» CPES- 


5% 


t 


D u -> K*~K~ 


0.380% 


* 


D u -> K*~K- 


0.220% 


* 


B- -> A^B U 


3.7 x 10" 4 


* 


B~ -> A*-B ,J 


6.1 x 10" 4 


* 


B" -» A-B* ,J 


3.6 x 10" 4 


* 


B~ K*~D*" 


7.2 x 10~ 4 


* 


B~ -> A"D U 


5.4 x 10" D 


t 


B~ -» A*-B ,J 


9.0 x 10~ D 


t 


B~ -> A~D* U 


5.3 x 10"° 


t 


B" -> K*~D*° 


1.06 x 10~° 


t 


B u ^ [KM ^ A-^+pO 


4.6 x 10 _b 


t 




5.5 x 10'° 


t 


B" - [A* u - A^+]B U 


6.0 x nr H 


t 


B" [A* u -> A-7T+]B* U 


7.0 x 10"° 


t 



TABLE V: Branching ratios for various B and D decays. Those from |T3 | are indicated by * while those which have been 
estimated have been indicated by f. 

III. CP ASYMMETRIES 

For a given decay of the form B~ -> K^~[T)^ -> F] where F is in general an exclusive or inclusive final state, 
let us denote 
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K+Tv- 




K + TV + 717T 


CPES- 


K*+K~ 


K*-K+ 


B- -> R-D 1 ' 


29 


53 


165 


1646 


70 


51 


B~ -> K*-D u 


28 


64 


219 


2636 


74 


44 


B~ -> i^IT u 


22 


47 


147 


1530 


66 


53 




40 


53 


244 


3590 


44 


45 



TABLE VI: The calculated branching ratio, d av — (d+d)/2 in units of 10 8 for the combination of the various 

decays in the left column with the D° decays in the top row using the parameters from Table llVl where the strong phases were 

chosen arbitrarily. As indicated, the results are averaged over each decay chain and its charge conjugate. 





K+TY- 


K* + -r- 


K + 1Y + 717T 


CPES- 




K*~K + 


B~ -> K-D u 


22 


12 


63 


905 


13 


9 


B- -> K*-D u 


21 


14 


83 


1450 


13 


8 


B- -> K-D* u 


17 


11 


56 


842 


12 


10 


B - _^ k*~D* u 


30 


12 


93 


1974 


8 


8 



TABLE VII: The branching ratio d av , in units of 10 8 , as in Table fVTl with the acceptance (not including R cu t) folded in. 



d = Br(B~ -» if (*)-[£)(*)" -» F] 
d = Br{B+ if(*)+[£)W° F] 



(3) 



where in the case of D*° it cascades down to D° before a D° — * F decay. In terms of the B and D branching ratios 
and the strong and weak phases, these quantities are given by 



d = acp + bcp + 2RF\facpbcF cos(Cs + Cf + t) 

d = acF + bcF + 2Rp \J acpbcF cos(£b + Cf — 7) (4) 

where a = Br{B~ -> K (*)-£)(*) ), b = Br(B~ -> K^-D^°), ( B is the strong phase difference of the B~ decay, 
is the strong phase difference for the D° decay and Rf is the coherence factor as defined in 0; in particular, for 
exclusive states Rf = 1. In the appendix, these expressions are generalized to the case where D°D° mixing is present. 
For each decay, let us define 



dav = (d + d)/2 

Acp = {d-d)/(d + d) (5) 

thus, Acp is the CP violating asymmetry. 
In Tabic EH 

we show the values of d av which result from the parameters of our sample calculation in units of 10 . 
In Table IVIIII we show the corresponding CP asymmetry Acp. In an ideal detector, the number of BB pairs required 
to observe the signal of CP violation (for definiteness, with a 3 a statistical error) is: 





R+TY- 


R* + tt- 


R + iv + mr 


CPES- 


R* + R~ 


R*-R + 


B~ -> K-D { > 


-58.2% 


-8.1% 


-17.9% 


-17.1% 


-17.1% 


14.3% 


B~ -> K*-D u 


-54.3% 


-63.3% 


-30.1% 


+11.6% 


-15.0% 


13.4% 


B~ -> R-D*" 


-76.9% 


-39.4% 


-30.2% 


-4.1% 


-12.1% 


8.1% 


B - K*~D W 


+58.4% 


-13.5% 


10.7% 


+15.8% 


+14.6% 


-7.2% 



TABLE VIII: The CP asymmetry Acp = (d — d)/(d + d) corresponding to the results in Table IVJ is given. 
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K+Tv- 




K + TV + 717T 


CPES- 


K*+K~ 


K*-K+ 


B- -> R-D 1 ' 


0.92 


25.88 


1.70 


0.19 


4.40 


8.63 




1.09 


0.35 


0.45 


0.25 


5.41 


11.39 




0.69 


1.23 


0.67 


3.50 


9.31 


25.88 


B~ -> K*-D w 


0.66 


9.32 


3.22 


0.10 


9.60 


38.58 



TABLE IX: The number N = 9/(d a vAc P ) of BB pairs needed to see a 3 — a signal of CP-violation in units of 10 s ; detection 
efficiencies, acceptances etc are not included in these numbers but are included in Table Ixl 









K + TV + 717T 


CPES- 


K*+K~ 


K—K + 


B~ -> K-D u 


6.05 


508 


22.4 


1.73 


122 


240 


B- -> K*-D u 


7.17 


7.95 


5.92 


2.27 


150 


306 


B- -> R-D* u 


4.53 


28.0 


8.82 


31.8 


259 


722 


B~ -» K*-D* u 


4.34 


211 


42.4 


0.9 


267 


1083 



TABLE X: The number N = 9/(d av Ac P ) of BB pairs needed to see a 3 — a signal of CP-violation in units of 10 8 taking into 
account the detection efficiencies, acceptances discussed in the text and assuming R cut = 0.2; see Eq. 



Tables IjlXIXIl indicate the number of BB pairs which would be required to observe a 3 — a signal of CP violation 
in these modes; Table lx) includes efficiencies and acceptance factors as given in Eq. From these tables it can be 
seen that with a (1 — 10) x 10 8 B's CP violation can be seen in several individual channels. In Table IXll we show the 
values of d av where we average over strong phases but keep 7 = 60°; likewise in Table IXlTl we show the r.m.s averaged 
values of |Ac7p|- Comparing these to the specific results in Table Ivnil obtained by assuming arbitrarily assigned 
values for strong phases, we see that the two approaches tend to give similar results. In Table IXlTTl we show the r.m.s 
asymmetries where we have also averaged over 7 where 7 ranges from to 2tt. 



IV. EXTRACTING 7 



We will now consider various strategies to determine 7 assuming we experimentally observe the results given in 
Tables <|VIIVIII|1 with the number of reconstructed events as discussed above. 

First, let us consider in isolation the case of B~ — * K~[D° — > K + tt~]. This rate, together with its charge conjugate 
gives us two distinct observables which are determined in terms of four unknown parameters: Ckd, Qk+tt-i b(KD) 
and 7. The two strong phases enter as the sum £ tot = C,kd + Ck+tt- so m effect there are only three parameters {(tot, 
b(KD), 7}. Still we cannot expect to reconstruct 7 but, as discussed in 0, this data gives a bound on sin 2 7. 

To illustrate this, in Fig. J5J) the thin solid line shows the minimum value of \ 2 as a function of 7. For each value 
of 7 we minimize with respect to the other parameters {(tot, b(KD)}. In the region of the graph where 7 > 45.6°, \ 2 
vanishes indicating that this is the lower bound within the first quadrant. Note that this and all similar graphs we will 
discuss are periodic with respect to 7 — * tt ± 7 so, for instance, in the second quadrant the bound is 7 < 180° — 45.6° 
etc.. 





K+TY- 


K* + -r- 


K + 1Y + 717T 


CPES- 


B~ -> K-D { > 


26 


38 


161 


1877 


B~ -> K*-D" 


43 


64 


268 


3095 


B- -> K~D W 


25 


38 


158 


1827 


B - k*~D w 


51 


76 


315 


3653 



TABLE XI: The value of d av is given using 7 = 60° and the branching ratios in Table IV! but with all of the strong phases taken 
at random. 
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K+Tv- 




K + n + nn 


CPES- 


Br -» K-D l> 


53.4% 


47.7% 


47.5% 


14.6% 


B~ -> K*~D U 


53.4% 


47.7% 


47.5% 


14.6% 


B~ -> R-D* u 


53.4% 


47.7% 


47.5% 


14.6% 


B~ -> K*-D* u 


36.8% 


33.0% 


23.0% 


10.4% 



TABLE XII: The r.m.s. average of |j4cp| corresponding to the results in Table fxTI 









K + n + nn 


CPES- 


Br -» K-D u 


44.3% 


39.3% 


39.1% 


12.0% 


B~ -> if*-D u 


44.3% 


39.3% 


39.1% 


12.0% 




44.3% 


39.3% 


39.1% 


12.0% 


S~ -> K*-D* u 


30.3% 


27.1% 


18.8% 


8.5% 



TABLE XIII: The results as in Table IXTTl averaged over all values of 7 (from to 27r). 



Clearly, in this case, the data is too meager to provide a useful bound. The 3cr bound (i.e. where \ 2 « 9) is only 
slightly above 0. 

We can also consider the bound on 7 obtained from the decay D° — > CPES— in isolation. There are more events 
of this type but the power this data to bound 7 is not much greater since Aqp is smaller (in general we expect the 
analyzing power of a particular mode to be ~ Aq P ). The minimum \ 2 m this case is shown with the dotted curve. 
Notice that taken in isolation it seems to be worse than even the single D° — > K + tt~ (CPNES) mode. 

Of course both of these two data sets depend on the common parameter b and if we have both sets of data together 
we obtain the results shown with thick solid curve which is an improvement on each of the data sets taken in isolation; 
in fact this thick solid curve gives a 3 a bound of 7 > 16°. As discussed in [f| since there the number of equations 
and observables is the same, there are ambiguous solutions which leads to the x 2 value being small over an extended 
range. 

To improve the situation, we can also use data from all four decays of the form B~ — > K^-*' . Note that 




gamma (deg) 



FIG. 2: The minimum value of x 2 is shown as a function of 7 for various combinations of data in the sample calculation. 
The thin solid line shows the result using just B~ —* K~\D° — » K + ty~] data. The dotted line shows the result using just 
the B~ -> is:"[D -> CPES-] data. The thick solid curve shows the result taking both B~ -> ^"[D -> K + n~] and 
B~ -> K~[T>° -> CPES-] data together. The dashed line shows the result using B' if(*)-[D ( * )0 -> K+n']. In the dash 
dotted curve, all four of the initial B~ decays where the D decays to the same two final states are considered. Thus the dashed 
dotted curve results from taking together data of the form B~ -> K ( * ] - [D ( * )0 K + n~] and B~ if(*)"[D ( * )0 CPES—]. 
The long dashed curve only includes data from two parent B~ decays, i.e. B~ — > K~D° as well as B~ — > K~ [D*° — > D°+7r°(7)] 
with either of the two D° decaying to K + ir~ as well as CPES—. The upper horizontal dashed line indicates the 3-sigma level 
determination of 7 with the luminosity required to give the results in Table [I] which corresponds to current B- factories while 
the lower horizontal dashed line corresponds to a luminosity 50 times greater which may be achieved at future high luminosity 
B factories. 
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FIG. 3: The minimum value of \ 2 using the four SCS modes in Eqn. J5J is shown with the solid line. The dashed dotted line 
shows the results with D° — * K* + K~ and D° — > K*~ K + together with the four parent B~ — ■> D ( * )0 _R" ( * )_ decays. 

each of these modes will have a different unknown value of b and C- In addition, the decay mode K*~D*° has three 
polarization amplitudes which we will take into account by introducing a coherence factor Rp into the fit since we 
are assuming that we are only observing the sum and we do not consider the additional information that could be 
determined from the angular distributions of the decays of the vectors as discussed in 0. If we consider the single 
decay D° — > K + tt~ we obtain the results shown by the dashed line which in this case gives a 3 — a bound on 7 of 
7 > 23°. The dot dash curve shows the result where we have both the D° — > K + ir~ and D° — > CPES— data. In 
this case we obtain a 3 — a determination of 7 (within the first quadrant) to be 60j^gg . Using the additional data 
improved the situation both by providing more statistics and because the different data sets have different spurious 
solutions leaving only the correct solution in common. 

For this dash-dot curve it is instructive to examine the number of observables versus the number of unknown free 
parameters. First of all, for D° — > K + tt~ there is the strong phase. For each of the four parent B~ decays there is a 
strong phase. In the case of B~ — > A*~D*° there is, in addition, a parameter R. Again, for each of the four parent 
decays there is the unknown branching ratio b = Br(B~ — > D*°) and finally the angle 7 giving a total of 11 

parameters. On the other hand, for each combination of B and D decays there are two observables, d and d giving a 
total of 16 so there is an overdetermination by 5 degrees of freedom. 

As another example, consider the case where only two the four combinations of B~ — > i^(*)~D(*)° are observed 
with D decay to K + tt~ and CPES—, then the system is still overdetermined. In Fig. [2] the long dashed line takes 
into account only the two B~ — > A"~D° and B" — > A~D*° and so has 6 unknown parameters determined by 8 
observables. Clearly having some overdetermination is helpful in obtaining a good determination of 7. 

It is important to contrast thick solid curve with the long dashed one, in Fig. [3 Recall both of them have 
D° — > A + 7r~, CPES—. However, in case of the thick solid curve the D° originate only from B~ — > A~D° whereas 
the long dashed curve is also getting the D° coming from D*° — > D° + 7T°(7). As a result whereas in the thick 
solid case there are 4 observables and 4 unknowns for the long-dashed case its 8 observables for 6 unknowns. That 
ends up making a significant difference as is evident from the figure; a lot more than one may naively expect just 
by doubling the number of D° or a factor of two in luminosity. Infact the lower horizontal line in Fig. |2 indicates 
(which corresponds to 3 — er determination at high luminosity) that the reduction in error on 7 is roughly a factor of 
5 (between the thick soilid and the long-dashed curves), i.e. a saving in effective luminosity of a factor of about 25. 

It is also instructive to compare dash-dot curve, which clearly has substantially more data, with the long-dash one. 
Notice that quality of determination of 7 by the two data sets is about the same. This suggests that once the number 
of observables is sufficiently large to overdetermine the parameters, further gains by including additional information 
only leads to modest gains. 

We can also determine 7 by using the four SCS modes in Eqn. (J2J. In Fig. [3] we show x 2 using these modes. If we 
again consider the four parent decays B~ — > DM /^ (*)~ then we obtain the results shown with the dashed dotted 
line which roughly gives 7 > 30° as a 3<r bound. 

Let us now consider additional sources of information which could constrain these results. In the next section we 
will discuss the impact of flavor tagging the D° by means of semileptonic decays. This is challenging but may be 
possible at B factories. Apart from B factory data, another source of additional information discussed in Q, is to 
use the charm factory to determine the strong phase differences Cf as well as Rf for inclusive modes. In Fig 0] the 
thick solid line shows the results using B~ — > K~D° followed by D° decay to A + 7r~ and CPES— as Fig EI The 
thin solid line shows the result if we suppose that we have the phase Ck+k- determined by a ip(3770) charm factory. 
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FIG. 4: The thick solid curve shows \ 2 as a function of 7 taking both B~ -> K~ [D° -»■ J^+tt"] and B~ -> X" [D° -> 
data together as in Fig. H The thin solid curve shows the result if Ck+tv- i s determined separately at a t/>(3770) charm 
factory. The dotted curve includes B~ -» if-[D° -> J^+tt"] B" -> K \D° -» ik*+7r-] B~ -> if - [D° -» ir+Tr" + nTr] 
J3~ — > A"~[D° — > CPES-] as well as the phase determination from a i/j(3770) charm factory. The dashed dotted curves 
considers all these decay modes together with the four parent decays B~ — > K*~ D*° as well as phase determination from a 
^(3770) charm factory. 
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FIG. 5: The results for B~ -» A"(*)~d(*)o f u owe d b y D° -» A'+tt"; CPES-; K+n'+mv, K'+n' where finite CP conserving 
backgrounds are included. The solid curve is for no background; the dashed curve is for signal/background=l, the dotted curve 
is for signal/background=i and the long dashed curve is for signal/background^^ 



Clearly this improves the situation by removing the ambiguities. Of course the situation can be improved still further 
by including several different D decay modes. In the dotted line, we show the result where we use ip(3770) data with 
the modes K + ir~ , K* + tt~ and CPES— and, in addition, the inclusive decay D° — > K + tt~ + mr. Note that as well 
as determining the strong phase difference for all the D° decays, the ?/>(3770) data also determines the value of Rf 
for F — K + n~ + mr. As before, we can also improve the situation by combining the data from all B~ — > Z)(*)° 
combinations as is shown by the dot dashed line. In this case, the determination of 7 (3cr) is 60 ± 10°. 

In the above, we have assumed that the there is no background. In general, of course, there should be a CP-even 
background which will tend to increase the error in 7. In Fig. 10 we show the x 2 as a function of 7 for all four modes 
including all the combinations of B~ — > Jf(*)~Z}(*)° with the solid line. The dashed line shows the results with a 
signal/background ratio of 1; the dotted line shows the result for signal/background is 4 while the long dashed curve 
shows the result for signal/background is 

The 3-cr errors in 7 in these cases is GOltil'o f° r no background and 60^26 1° f° r signal/background=l. 

Recall that the color suppressed branching ratios for B~ — > K~ D° and B~ — ► K*~D° are likely to be very difficult 
to measure experimentally (see however section EJ. The above analysis yjul is therefore designed to yield both these 
suppressed Br's as well as 7 with the input of each data set. Fig. serves to illustrate how well this works out 
using data for all the combinations of B~ -> lf(*)-D<» followed by D° -> K+n-; CPES-; K + tt~ + mr, K*+n'. 
The dependence on Br(B~ — > K~D°) is indicated by the solid line while the dependence on Br(B~ — > K*~D°) 
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FIG. 6: The minimum value of \ 2 is shown as a function of Br(B~ — > K~ D°) and Br(B~ — > K*~ D°) using data for 
all the combinations of B~ -> A' ( * )_ D W ° followed by D° -> K + tt~- CPES-; K + ty~ + mv, K* + -k~ . The dependence on 
Br(B~ — > K~ D°) is indicated by the solid line while the dependence on Br(B~ — > K*~ D°) is shown by the dashed curve. 



is shown by the dashed curve. The 3-sigma determination of Br(B — > _ftT D°) is thus » ±1.8 x 10 6 and for 
Br(B~ -> K*~D°) it is » ±2.25 x 10" 6 . 

V. FLAVOR TAGGING D° MESONS VIA SEMILEPTONIC DECAYS 

Clearly a very useful constraint on the system of equations Eqn. Q would be a direct determination of the branching 
ratio, b = Br(B~ — > K~D°). This can only be done if one observes the I) decay to a flavor specific semi-leptonic 
decay 0. Unfortunately this is subject to a large background from semi-leptonic decays of the parent B~ since both 
semi-leptonic decays result in negatively charged leptons. 

The kinematics of the two kinds of processes are, however, quite different so it may eventually be possible to 
determine b. As a case in point let us examine the kinematics involved if one considers B — » K~D° followed by 
D° -> £-9 e K+. 

First of all, X^= e u Br(D° — > (,~v?K) = 6.9%, hence the combined branching ratio using the estimate in Table IVI 

is: 

Br(B~ -> K-[D° -> fv t K + )] = 3.7 x 1CT 7 (7) 

The final state for this signal is thus vi~K + K~ and it will be subject to the following backgrounds (including 1 = 
e and pi) : 

(a) The semi-leptonic decay B~ — > D£~[D° — > K + K~\. This has a combined branching ratio of 1.7 x 10 -4 . 

(b) The semi-leptonic decay B~~ — » Di~[D° — > 7r + i^ _ ] where the 7r + is misidentified as a The combined 
branching ratio is 1.6 x 10~ 3 . Of course this must be multiplied by the rate of mis-identification for a given 
state. 

(c) B~ -> p£~ + [X c -> if+if- + X], with a branching ratio » (10^ 2 - 10" 3 ). 

(d) £? - — > D£~K + K~; the branching ratio is crudely estimated to be around 3 x 10~ 5 . 

(e) Backgrounds form continuum events 

Against these backgrounds one can apply the following kinematic cuts. 

(1) For the signal the energy of the K~ in the rest frame of the B~ is fixed to be E K - = (m 2 B + m 2 K — m 2 D )/ (2ms). 

(2) The missing neutrino leads to the constraint \ps — Vk~ ~ Pi ~ Pk+ I = 0. 
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FIG. 7: x 2 as a function of 7 for B — > K [D° — + CPES— } . The dotted curve is the same as in Fig. |5J the long dashed curve 
assumes that 6 has a 10% 1-a error, the thick line assumes a 100% error and the thin solid line assumes a 200% error. 

(3) For backgrounds of the type (a), the invariant mass of the K + K~ system will be mn; indeed, to eliminate the 
related backgrounds from D° — > K + K~ + X one may use the cut \pk+ + Pk~ I > m D- 

(4) Signal events will have three distinct vertices while continuum background will have only 1. 

Background (a) and (b) are ~ 10 3 x signal. Cut (3) would remove these entirely except for momentum resolution 
but the additional cuts (1) and (2) which are satisfied by the signal may be sufficient to control this large background. 
Backgrounds (d) is probably O(10 2 ) times the signal; (c) is probably even a lot bigger, (d) in particular would pass 
cuts (2) and (3) and would have to be reduced by cuts (1) and (4). Likewise background (c) would have to be reduced 
by the monoenergetic kaon cut (1). Background (e) could, in principle, also be large and it is not clear whether cuts 
(1) and (4) together with the standard cuts against the continuum can be sufficient to eliminate it. 

In spite of the difficulty in determining b, in the case where the D° decays to a CP eigenstate, a relatively weak 
bound, on b, can be helpful in a determination of 7. 

In Fig. [3 we show the x 2 curve obtained just from the GLW modes, B~ — » K~[D a — ► CPES-], as in Fig. J2J with 
the dotted line. The dashed line, the thick solid line and the thin solid line correspond to the result obtained assuming 
a 10%, 100% and 200% 1-a Gaussian error in the determination of b respectively. These lead to 3-<r lower bounds on 
7 of 25.2°, 16.8° and 12.6° with the prospect of improvement at higher statistics. The reason for this is that the CP 
violation in this case is relatively small ~ 15% so that the solutions with small 7 correspond to a situation where b 
is large but Cb ~ 180° so there is near cancellation between the two channels. This is quite different from the actual 
situation where b << a so a modest bound on b improves the situation greatly. 

If an experimental value for b is not available, a theoretical estimate can serve at the expense of the result becoming 
model dependent. From the above, such an estimate need not be very precise to be of some utility. On the other 
hand, the decay B~ — > K~D° is color suppressed making reliable theoretical predictions more difficult. 

In Fig. [8] we show the analogous calculation where we have combined the data from all the combinations of B~ — > 
#(*)-[£)(*)o CPES-] with a determination of b to 10% (thick solid line), 200% (dotted line) and unconstrained 
(dashed line). These give the corresponding 3-sigma lower bound on 7 of 40.2°, 16.2° and 3° respectively and upper 
bounds of 76.8°, 78.0° and 78.6°, respectively. 

VI. IMPACT OF DD MIXING 

As has been discussed previously in the literature 0, , the effects of DD mixing could be significant on some 
of the combined branching ratios. This is particularly true for final states such as K + ir~ because the favored decay 
B~ — > K~D° could be followed by the favored D° — ► A + 7r~ if the D° were to oscillate to a D°. If the probability 
of oscillation is 0(1%) then this channel might be comparable to the direct DCS rate for B~ — > K~[D° — ► A + 7r~] 
which is assumed to be the only channel in the absence of oscillation. 

As discussed in the appendix, mixing is controlled by 4 parameters, x, y, A and <p. If CP is conserved in D decays, 
which is an excellent approximation in the SM, then A — <fi — 0. The standard model predicts that the parameter 
controlling the mass difference, x ~ O(10~ 4 ). The Standard Model prediction for the parameter relevant for the 
life time difference, y, is less certain since it could be enhanced by final state hadronic interactions. It has been 
suggested j2^ that y could be as large as 10~ 2 . 
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FIG. 9: The solid line shows the fit to the data from K + tv~ , K + ir~ + nir, CPES— and K* + ix~ with no D mixing, the dashed 
line shows the fit, assuming y = to the data in the case y = 0.01 the dotted line shows the fit, assuming y = to the data 
in the case y — 0.1 while the dash dotted line shows fit where a non-zero value of y is taken into account where the data is 
generated with y — 0.1. 



In Fig. |21 the solid line is a fit to the data for K + tc~, K + ir~ + mr, CPES— and K* + ir~ final states as discussed 
before. The dotted line shows the result if y — 0.01 but the result is fit assuming there is no mixing. The minimum 
value of x 2 has been shifted down to 58°. The dashed line shows the same calculation for y — 0.05 and the minimum 
value of x 2 is shifted down to 51°. The dashed dotted line shows the x 2 f° r V — 0-05 where the fit is done taking into 
account a non-zero value of y. In this case the fit value of 7 is 60l^g' 5 . 

Of course data from a ^(3770) charm factory should be able to improve significantly on the constraints on the 
D°D° mixing parameters and thereby it could immensely aid in the accurate determination of 7 from the B — > KD° 
method. 



VII. SELF TAGGING NEUTRAL B DECAYS 



It is also useful to add information obtainable from the decay of neutral i?'s. Since B° undergoes oscillation there 
are two strategies which can be considered: 

• The B° decays which are flavor non-specific such as B° — > KsD° or B° — > K*°D° where K*° — * Kgir . In this 
case the oscillation of the B° plays a role. 

• Self-tagging 0,0] decays which are flavor specific such as B° — * K*°D° where K*° — > K ~n + . In this case the 
oscillation of the B° does not play a role. 

In the flavor non-specific case, the oscillation phase (3 in the Standard Model enters in such a way that the main 
dependence is on 8 = 2/3 + 7 = /3 — a + ir. This has been discussed in [T^. IT3I llH HH . 
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FIG. 10: The solid curve is generated using charged B~ data and is the same as the solid curve in Fig. The dashed curve 
uses only neutral data with the two parent decays B° — > [K*° — > _K" _ 7r°]D'*- )0 ] and subsequent D decay to K + tt~, K* + tt~, 
K+ n + nn, CPES—. The dotted curve combines the data set used in the solid curve (charged B) with the data set used in 
the dashed curve (neutral B). 

Here we will discuss the former case where the oscillation plays no role and so only direct CP violation via flavor 
specific decays is used in a determination of 7; for instance, via the decay B° — > [K*° — > _ftT~7r + ]D°. The main 
difference with the charged (B^) decay is that the b — > c channel is color suppressed and thus will be about 9 times 
smaller than the charged case. Thus we will estimate the branching ratios as shown in Table Ivl The b — > u channel 
proceeds via the same diagram in the charged and neutral cases hence the branching ratios in these cases should be 
about the same as shown in Table Ivl 

Because of the suppression of the b — ► c channel, if D° decays to a DCS mode the amplitudes no longer are matched 
and so the CP violation will be somewhat smaller than in the charged case. Conversely the CP violation in modes 
with D decays to CPES will now be larger. 

In Fig. [23 we illustrate the use of such modes. The solid line, for reference is the same as in Fig. 03 including data 
from all combinations of B~ -> K^~T)^° with subsequent D decay to K+n~, K* + n~, K+tt- + nn, CPES-. The 
dashed curve considers the neutral B decays B° — ► [K*° — > K~n°] with subsequent D decay to the same four final 
states. Because of the extra color suppression, it is not surprising that it falls below the solid line. Of course we 
can combine the two data sets since the strong phases of the D decays should be common to each in which case one 
obtains the dotted curve improving the situation slightly over the solid curve. 

VIII. MULTIBODY B DECAYS 

Two body decay modes of the B to D° suffer from the problem that at least one of the channels is color suppressed. 
A method to circumvent this has been suggested by Aleksan, Petersen and Soffer (APS) 0] that decays of the form 
B~ — > K~n°D° may be used to extract 7. On the quark level the processes are the same as the two body B decays 
considered above but in this manifestation there are two important differences: 

1. Both the b — > c and b — + u channels are color allowed. 

2. The amplitude is a function of the 3-body phase space. 

The fact that the b — > u channel is color allowed is important since it can improve the statistics, particularly if D 
decays to CPES. To take full advantage of this, however, we must overcome the fact that this is an inclusive state 
spread over phase space (i.e. the Dalitz Plot). 

The APS method 0] is to extract from the data a number of different amplitudes and their relative phases by 
fitting to the Dalitz plot using a model with resonance and continuum terms and comparing the D° decay to CPES 
and flavor specific modes. 

Here we will consider a complimentary approach where we integrate over the Dalitz plot using the formalism of 
Eqn. (|13|) in the Appendix where there is an additional coherence parameter associated with the B decay. In this 
approach, though, several different three or more body decays may be taken together which will tend to increase the 
statistics. 
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FIG. 11: The solid curve is the same as in Fig.Q)]while the dashed one is the \ 2 f° r 7 determined by the three body B decay 
modes D°_R-7r°, D* _R"7r° and D°K*- ~7T° using the assumptions discussed in the text. 



First, let us compare the number of degrees of freedom versus the number of observables in such a case. Suppose 
that we consider N inclusive decays of the form B — ► T)°K + Xi (i = 1 . . . N) where the D° subsequently decays into 
CPES-, Me exclusive modes and Mj inclusive modes. For each B decay there are 3 parameters, b, £ and Rf- For 
each inclusive D mode there is one phase while for each inclusive D mode there is a phase and a coherence factor. In 
addition, of course, 7 is also unknown giving 3 N + Me + 2Mj + 1 parameters. On the other hand, for each combination 
of B decay modes with each D decay mode there are two observables, d and d giving a total of 2N(M E + Mj + 1) 
observables. Thus, 7 may be determined if 



2N(M E + Mi + 1) > 37V + M E + 2Mj + 1 (8) 

where for the case of equality applied there would likely be ambiguities. If we rearrange the above equation assuming 
that M E + Mi > we obtain: 



n> Me±2M l±1 

~ 2M E + 2Mj - 1 w 

From this form we see that the lower bound on TV is 3 only in the case of Me = and Mi = 1. In other cases the 
lower bound is either 1 or 2 (taking into account that N is an integer). Thus at least 1, 2 or 3 different B decays 
must be observed depending (according to eqn. JJjJ) on how many different D decays are considered assuming, in all 
cases, that we observe the CPES- decays of the D°. 

As an illustration of this let us assume that the three inclusive modes B~~ — > D ^ - ^ , B~ — * D*°X~7r° and 
B~ -> D°K*-n° are observed with D° -> CPES—, K+tt- and K*+t:- so that N = 3, M E = 2 and Mj = 0. Clearly 
Eqn. (JHJ is satisfied since there are 18 observables for 12 unknowns. To carry out our calculation we will assume that 
a(Z)(*)°Jf(*)-7r ) = a(D^°K W~) and b(DW Q K*-ir ) = 9b{D^°K^-) since it is not color suppressed; also we take 
R = 0.5 in each case. The results are shown with the dashed line in Fig. 1111 



IX. PROSPECTS FOR 7 AT A SUPER-B FACTORY 



Improvements at the existing asymmetric B-factories at KEK and SLAC are expected to allow luminosity increase 
to around 10 35 cm~ 2 sec~ 1 . Further increase in luminosity to 10 36 cm -2 sec -1 and beyond will most likely require 
a new machine 27]. The asymmetric B-factories at KEK and SLAC have made a rather accurate determination of 
sin(2/3) and established that the CKM- phase is the dominant contributor to the observed CP asymmetry in B — > ipK s , 
which has to be regarded an important milestone in our understanding of CP violation phenomena. There is now 
considerable interest at the construction of an asymmetric "Super-B Factory" (SBF) with luminosity 10 36 cm~ 2 sec~ 1 
or more pRll^ . Precision determinations of all the angles of the unitarity triangle, roughly in the range of the intrinsic 
theory error for each of the three angles, and the search for beyond the SM source(s) of CP violation, which are bound 
to exist, constitute important motivations for such an effort. 
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In particular, 7 can be determined very cleanly from direct CP violation via the generic B — > K(*)D°(*) processes. 
There are multitude of available strategies and modes; several of them we dealt in this paper. Using Fig[3]as a guide 
we can anticipate possible determination of 7 at a Super-B factory with a 3-er error of a few (as 2) degrees. 

Admittedly backgrounds could make things difficult in an actual experimental setup, as Fig indicates; from the 
figure we see that the accuracy in determination of 7 from these methods may, due to backgrounds, get degraded 
to w 7° at the Super-B Factory. It should be noted though that strategies used in Fig [5] and Fig [5] are far from 
exhaustive. For example, one promising technique which we had suggested earlier 00 is Dalitz plot study of 3-body 
decays of D°. Since some experimental data had existed at that time for D° — > K + ir~ir° from Fermilab experiment 
E687 123, we had used that specific mode to illustrate the method. More recently CLEO collaboration |3 1) has also 
studied D° — * K s tt + ti~ which is another interesting candidate where the analysis could be readily applied to 
extract 7. Indeed preliminary studies of this mode by the BELLE collaoboration |32j are quite encouraging. 

In addition to Dalitz plot studies of D° decays, charm factory running on -0(3770) could be very helpful in deter- 
mination of 7 from B — > KD approaches as Fig Q] illustrates; see also section IVIIII and 7]. 

X. SUMMARY AND OUTLOOK 

In this paper we have extended our previous studies 0,0,0 on extraction of 7 using direct CP violation in B — > KD 
processes. In principle, these methods are theoretically very clean. The irreducible theory error originating from 
higher-order weak interactions is O(10 -3 ) 26], i.e. in all likelihood even smaller than the theory error in deducing the 
angle (3 using time dependent CP asymmetry in B — > ipK s . However, 7 determination from B — > KD is much harder 
than /3 from B — > ipK a . 

This study strongly suggests that the demands on machine luminosity can be significantly alleviated if a combination 
of strategies is used. One interesting handle that we examined here which looks rather promising is to include D*° 
from B -> K^D*°. The formalism for the use of D° decays is identical to B -> K*T>° after D*° -> vr°(7) + D°. 

Similarly including K* (via e.g. B~ — ► K*~D°) along with B~ — > K~D° is helpful. Also it of course helps a 
great deal to use both CP eigenstatesQ along with CP non-eigenstates of D°, whether they be doubly Cabibbo 
suppressed00 or singly Cabibbo suppressed \v?\. 

While flavor tagging of the D° and D° is likely to be a challenge, we emphasized that kinematics and topology 
of the signal events is quite distinct from that of many of the backgrounds so that semi-leptonic tags may have a 
chance at least for a relatively imprecise determination of the needed color-suppressed branching ratio. If such a 
determination could be made, say at super-B factory with an error on the branching ratio of around 10-20%, it could 
become helpful in 7 determination. 

Many of the modes relevant for gamma determination should exhibit large direct CP-asymmetries; this is especially 
so in the case of doubly Cabibbo suppressed modes of D°. However, as Tables (|VIIIIXIIIXIII(1 illustrate CP asymme- 
tries in many other modes are also appreciable. Indeed, at least in the initial stages, it may even be useful to target 
the search for such large asymmetries. 

It is difficult to overemphasize the important role that a charm factory running at ^(3770) can play in the extraction 
of 7. Charm factory can of course help by improving on the mixing parameters of the D° — D complex. It can 
also help in determination of the doubly Cabibbo suppressed branching ratios as well as in the needed strong phases. 
Quantum entanglement in tp — > D°D° decays can also be exploited to determine the coherence factor and the strong 
phase for inclusive D° decays 0. 

Note also that all of these studies of charm physics at the charm factory not only have important application to 
extraction of 7 from direct CP studies in B — > KD processes, precisely the same information on D-decays and D- 
mixing can also be used in time dependent CP studies in B° — > K°D°(*); in that case one gets the linear combination 
of unitarity angles 8 = (2/3 + 7) = (0~ ot + n) and furthermore that method also gives /30,Q providing an additional 
test of the CKM-paradigm. 

While B-factories with about 10 9 B-pairs are likely to be able to make appreciable progress in determination of 7, 
super-B factory with luminosity > 10 10 B-pairs will be needed to extract 7 with an accuracy roughly commensurate 
with the intrinsic theory error that these methods allow. This in itself should constitute an important goal of B-physics 
in general and super-B factory in particular. 

XI. ACKNOWLEDGMENTS 

We would like to thank Tom Browder, Riccardo Faccini, Tim Gershon, Masashi Hazumi, Gerald Raven, Abi 
Soffcr, Sheldon Stone and Sanjay Swain. This research was supported by Contract Nos. DE-FG02-94ER40817 and 
DE-AC02-98CH10886. 



18 



Appendix: Formalism 

Let us now present a summary of the formalism for determining the rates for B~ — * [£)(*)° — > F] and related 

processes. Here we will discuss the rates in four situations: 

1. F is an exclusive state; without D°D° oscillations. 

2. F is an inclusive set of states; without D°D° oscillations. 

3. F is an exclusive state where D°D° oscillations are considered 

4. F is an inclusive state where D°D° oscillations are considered. 

For this purpose exclusive state refers to a state which is governed by a single quantum amplitude, for instance 
D — > K + ir~. In contrast an inclusive state is composed of a decay where there are a number of amplitudes 
contributing incoherently. Some examples are inclusive states made of states with different particle content, for 
instance K + + mr or multibody final states integrated over phase space, for instance K + tt + tt° integrated over the 
Dalitz plot. Note that each point in the K + tt~tt° Dalitz plot is an exclusive state. 

First let us take / to be an exclusive state. Let us denote by d the branching ratio for B~ — ► K^*' [Z)(*) — > /] where 
D° means a general mixture of D° and D° and J the charge conjugate branching ratio for B + — > K^ + [D^° — > /]. 

For this decay let a be the branching ratio for B~ — > K~D° and b be the branching ratio for B~ — > K~D°. 
Likewise let us denote by c the branching ratio for D° — ► / and c the branching ratio for _D° — * /. Assuming that D° 
mixing is ncglig ibleHH, 



d = acf + bef + 2a/ abcfCf cos(( B + Cf +l) 

d = acf + bef + 2yJabcfCf cos(Cs + Cf ~ l) (10) 

where Cb is the strong phase difference between B~ — > K(*^~ £}(*)° and B~ — > K(*)~ f)(*) and the strong phase 
difference between £>° — > / and £)° — > /. 

Let us now consider the generalization to the case where F is an inclusive state, F = {/i, . . . , /„} and denote: 

Jfce«* = ^-^g^ 6 ^ (11) 

where c F = J^. cj 4 and similarly . For the decay to an inclusive final state F, then, we can generalize eqn. I|1U|) 
to: 



c? = acF + &Cf 

+2i?^\/ abcFCF cos(Cb + C-F + 7) 
d = acf + 6c> 

+2R F ^/abc F c F cos(Cs + Cf - 7) (12) 

Note that in the case where F is an exclusive state R F = 1. 

Another, related generalization is the case where the initial £>~ decays to an inclusive state. For instance this is 
the case in B~ — * K*~D*° because there are three helicity amplitudes. Of course the amplitudes may be separated 
through the consideration of angular distributions in the vector decays. Another instance is a three body decay such 
as B~ -> KnD°. 

In general, we can consider B~ — > G where G is an inclusive set of states. Thus, each gr, contains one (generic) D° 
meson. We will use g[ to denote, specifically, the version of the state where the neutral D is in a \D° > state and 
g'l to be the case where the neutral D is in a \D° > state. Let us denote the branching ratios B~ — * g\ by a gi and 
B~ — > g" by 6 9i . One could, in principle, determine a gi and b gi by observing the semileptonic decay of the D . 

Thus, in analogy to Eqn. [fTT|l let us define: 



R G e lCG = 



V a 9. b 9 



(13) 
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If the neutral D in G subsequently decays to an inclusive state F, then the combined branching ratio is: 



d = <igcf + baCF 



+2R F R G ^a G b G c F c F cos(Cg + Cf + 7) 
d — a G c F + b G c F 

+2R F R G \J a G b G c F c F cos(Cg + Cf - 7) (14) 
Let us now consider the case where D°D° mixing is present. 

First, let us consider the fully general case where, following the usual formalism 0,113, the eigenstates of neutral 
D meson are: 



|Z> X ) = p\D°)+ q \D ) 
\D 2 ) = p\D°)-q\D ) 

with corresponding masses m\, mi and widths Y\, Y-i. These are characterized by the parameters: 



(15) 



771 2 — 777 1 



y 



2r 



(16) 



where 2r = T 1 +T 2 . 

If p/q 7^ 1 then CP is violated in the mixing. If this is indeed the case, one would expect that there would also be 
CP violating phases in the decay. In the following we will assume that D decays are CP conserving and then indicate 
how to generalize to the case of CP violation in the final state. In this case the decay amplitudes of D° only have 
strong phases and we denote: 



V - = exp A+i * 

q 

Thus, generalizing Eqn. (|14(l to the case where there is D° mixing and integrating over time: 

, Q + P\ , . - 

d = — - — a G c F + b G c F 

+2R F R G y / a G c F b G c F cos(?7 G + 7 + rj F ) 



(17) 



Q-P 



OA v a . 

e a G c F + e b G c F 



-2Ai 



-2R F R G \J a G c F b G c F cos(rj G + 7 - r/p) 



-yQ c F \/ 'a G b G R G e A cos(t7 G + 7 + 

+c F V a G b G R G e~ A cos(r] G + 7 - 4>) 
+a G y/c F CFe~ Rf cos(i]f + 4>) 

+b G ^CFCFe Rf cos(r]F — <5 



+xP 



c F \J a G b G R G e A sin(?7 G + 7 + , 



+c F \/ a Gb G R G e A sin(?7 G +7- 
+a G y/cFCFe~ A Rf sm(r) F + 0) 

+b G \J CFCFe A Rf sm(rjF — (, 



(7) 



Q + P 



a G c F + b G c F 



+2R F R G y / a G c F b G c F cos(ry G - 7 + rj F ) 



Q-P 



e agcj? + e b G c F 



-2R F R G \J a G c F b G c F cos(r) G - 7 - r) F ) 



-yQ 



c F \/ a G b G R G e A cos(rj, 



'G - 7 ' 



+cf y/a G b G R G e A cos(r] G - 7 + 
+a G ^/c F c F e A R F cos(rj F — 0) 

+b G i/c F E F e~ A R F cos(rj F + <, 



where 



Note that in the limit that x, y, A, 



-xP 



cf\/ a G b G R G e A sin(r] G 



-c F y/a G b G R G e A sin(7y G - 7 
-a G \J c F c F e A R F sin(?7F — 4>) 



b G \J c F c F e R F siri(i] F + < 



P 



1 



Q 



1 



1 + x 2 ' 1 — J/ 2 

these equations reduce to Eqn. Ijl4(l . 
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